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Abstract 

We present a family of models for term structure dynamics in an 
attempt to describe several statistical features observed in empirical 
studies of forward rate curves by decomposing the deformations of the 
term structure into the variations of the short rate, the long rate and 
the fluctuations of the curve around its average shape. This fluctua- 
tion is then described as a solution of a stochastic evolution equation 
in an infinite dimensional space. In the case where deformations are 
local in maturity, this equation reduces to a stochastic PDE, of which 
we give the simplest example. We discuss the properties of the so- 
lutions and show that they capture in a parsimonious manner the 
essential features of yield curve dynamics: imperfect correlation be- 
tween maturities, mean reversion of interest rates and the structure 
of principal components of term structure deformations. Finally, we 
discuss calibration issues and show that the model parameters have a 
natural interpretation in terms of empirically observed quantities. 



Nous proposons une description de la deformation de la structure 
par terme des taux d'interet en termes d'une equation aux derivees 
partielles stochastique. En etudiant en detail le cas lineaire, nous 
montrons qu'une telle formulation rend compte des observations em- 
piriques sur la dynamique de la courbe des taux forwards avec un 
nombre restreint de parametres. Nous discutons en detail le role des 
hypotheses mathematiques et leur impact sur la dynamique des taux 
d'interet. 
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1 Introduction 

1.1 Motivations for term structure modeling 

There are two very different motivations for term structure modeling. Tlie 
first is concerned witfi tlie pricing of interest rate derivative securities. In 
tliis context, whicli fias been tfie principal motivation behind term structure 
models in the mathematical finance literature |piR| , |HJM| , |Vasicek|| , the main 



concern has been the development of "coherent" -in the sense of arbitrage- 
free - pricing criteria for securities whose payoffs depend on movements of 
interest rates. 

The second motivation, which could be labeled as "econometric", is the 
statistical decription of the movements of real interest rates. In contrast to 
the preceding approach where the emphasis is on cross- sectional coherence of 
prices given by the model, here the emphasis is on describing and reproducing 
as closely as possible the time evolution of interest rates from a statistical 
point of view. Such an approach is useful if one is interested in simulating 
scenarios, calculating Value-at-Risk of fixed-income positions but also from 
a theoretical point of view, to gain a better understanding of interest rate 
fluctations and their relations to other economic variables.. 

From a mathematical point of view, the first approach corresponds to 
modeling the dynamics of interest rates under a risk-neutral (or risk-adjusted) 
measure, while the second approach corresponds to the modeling of real- world 
term structure dynamics. 

Given the complexity of the behavior of the yield curve, there is a conflict 
between these two criteria which turn out to be difficult to satisfy in the same 
model. Most of the existing work in the mathematical finance literature on 
interest rate modeling has adopted the first approach, namely the description 
of yield curve dynamics under a "risk- adjusted" measure. The reason for this 
trend is not difficult to understand: contrarily to a stock market, in a bond 
market the future values of many securities - namely the zero-coupon bonds - 
are known with certainty. The observability of the contemporaneous prices of 
these bonds then makes it possible to cahbrate a model for the risk-adjusted 
dynamics of interest rates directly to the observed bond prices. However, as 



correctly pointed out by various authors (see e.g. fPliska 1997|| , pp. 201-204), 
once the risk-adjusted dynamics has been calibrated it is not obvious that 
such a model will tell us anything useful about the real-world dynamics of 



interest rates. 

In fact, it turns to be difficult to account for a series of empirical facts ob- 



served in studies on the term structure of interest rates llBouchaud et al., 1997 



■) 



in the framework of standard term structure models of the first type. On one 
hand, it seems that the constraints implied by the absence of arbitrage in 
these models are so strong that the latter is obtained at the detriment of a cor- 



rect representation of the dynamics of the yield curve jBouchaud et a/., 1998|, 
El Karoui et a/.,1996|| . On the other hand, some stylized empirical facts 
about term structure deformations seem to have no theoretical counterpart 
in classical arbitrage-free models. 

1.2 Continuous models for discrete observations 

Empirical term structure data consist of time series of interest rates of various 
maturities; for example, in the Eurodollar market interest rates of about 40 
different maturities can be obtained on a daily basis. Instead of modeling the 
data as a sequence of 40-component random vectors, theoretical models tend 
to represent the interest rate curve as a function of a continuous maturity 
variable 6 evolving with a continuous time parameter t. These models, which 
offer greater analytical tractability, correspond to the ideal case where fixed- 
income instruments for continuum of maturities are traded in a continuous- 
time market. Continuous models for interest rate dynamics can be divided 
into two classes. 

The ffist class starts by postulating a stochastic model for the short term 
interest rate: typically, the short rate is considered to follow a stochastic dif- 
ferential equation. The model then specifies a rule for constructing the yield 
curve from the short rate by an additional set of variables representing the 
market price(s) of risk. This class includes the Vasicek model ||Vasicek|| , the 



Cox-Ingersoll-Ross model ||CIR|| and other variants of these models. Short- 



rate based models present the advantage of analytical tractability. However, 
they reduce the dynamics of the whole yield curve to the movement of one 
it its endpoints, which makes these models rather rigid. A more serious 
objection to short-rate based models is that empirical studies tend to re- 
ject the hypothesis that the short rate is Markovian, let alone a diffusion 
|Ait-Sahalia 1997|| . 



The second class of models, initiated by Ho & Lee [[Ho fc Leej and ex- 



tended by Heath, Jarrow & Morton |[HJM|| , take the initial term structure 



and the forward rate volatilities as inputs, model the evolution of forward 
rates as a family of scalar diffusion process indexed by the maturity date T: 

dfK3Mit,T) = a{t,T)dt + ait,T)dWt (1) 

The forward rate curve is therefore a continuum of diffusion processes driven 
by a finite number of noise sources; it is therefore intuitively clear that if no 
constraint is imposed on the drift and diffusion coefficients such a model will 



present obvious arbitrage opportunities. It was shown in ||HJM|| that, when 
the noise source is a finite dimensional Wiener process, then the absence of 
an arbitrage strategy involving bonds of all maturities imposes the following 
relation between the drift and the diffusion process of the forward rates: 

a{t,T) = (j{t,T) f a{t,u)du + a{t,T)-f{t) (2) 



where 7(t) is some predictable process independent of the time to maturity 
6 = T — t. Underlying this result are the hypotheses that: 

- the number of sources of randomness is finite (in practice, much smaller 
than the number of maturities) 

- arbitrage opportunities involve any number of maturities. 

The first hypothesis is a modeling choice which means that for repre- 
senting the fiuctuations of the yield curve, a random vector with around 40 
components, we use a model where the number of components A^ is infinite 
(in fact a continuum) but the number k of sources of randomness is kept 
finite. The well-foundedness of this choice is not obvious a priori and jus- 
tified by observing that Principal Component Analysis (see below) of yield 
curve movements reveal that around 95% of the variance of bond yields is 
explained by 3 factors. As we shall see below, such an observation does not 
imply that properties other than the covariance structure will be explained 
by the same small number of factors. Moreover, the relation between drift 
and volatility (Eq.§ ) which is fundamental in HJM models is verified if the 
sources of randomness defining yield curve dynamics are correctly specified. 
What if one neglects a source of randomness when specifiying the model? 
This amounts to "projecting" the model on the space spanned by the first 
few factors or integrating the drift and volatility functions in Eq.(|I|) with 
respect to the left-out variable. In this case relation (Q) has no reason to 
hold anymore between the modified drift and volatility. In other words, the 

6 



HJM relation (|I]) is not robust to misspecification of the number of sources 
of randomness. 

An alternative approach, which we propose here, is to consider the limit 
where k is also infinite, thus restoring more flexibility to the forward rate 
curve and at the same time reducing the possibility of easy arbitrage [|. At 
the same time, we shall take into account the relative segmentation of the 
fixed income market into maturity-specific markets to introduce the notion 
of local deformation and examine its consequences, as explained below. 

1.3 The forward rate curve as an infinite-dimensional 
process 

In many mathematical models of the yield curve, the term structure of inter- 
est rates is often parametrized through the forward rate curve. Let B{t,T) 
be the price at time T of a (default-free) zero-coupon bond with maturity 
date T. The instantaneous forward rates /(t, 6) are related to bond prices 
by 

Vte [0,T],5(t,T) = exp(- / f{t,u) du) (3) 

JO 

The forward rate term structure at time t is therefore a function of a con- 
tinuous variable 9, representing the time to maturity. As such, the forward 
rate curve /(t, .) naturally lives in an infinite-dimensional space of functions 
and, as remarked by several authors [ Dufhe fc Kan , Kennedy 1994 1 there is 



no reason a priori to believe that its random movements will be restricted 
to a finite-dimensional subspace or manifold. There are two main reasons 
why the existing theory has almost exclusively focused on models with a fi- 
nite number of state variables: the first reason is analytical tractability, the 
second is that principal component analysis of term structure movements 
seems to suggest that a few factors are sufficient to explain the covariance 
structure of interest rates [[Bouchaud et ai, 1997] , [Litterman fc Scheinkmaii 



As pointed out above, making a large (formally infinite) number of forward 
rates dependent on a small number of factors restricts considerable the vari- 
ous configurations of term structure movements, creating a confiict between 
tractability and a faithful representation of empirical observations. 



^On arbitrage in infinite factor models, see [ Douady[ 



Note however that, as pointed out by |[Musiela|| , even with a finite-dimensional 
noise source in a HJM approach, one cannot ignore the infinite-dimensional 
character of the forward curve process. This means that, in most such mod- 
els, the smallest family of forward rates which is Markovian will be an infinite 
family. However, HJM models endow this infinite family of forward rates with 
only finitely many degrees of freedom, placing considerable restriction on the 
type of dynamics it can follow. The difference in the present approach is not 
only to consider the yield curve as an infinite-dimensional process but also 
as a process with infinitely many degrees of freedom. 



1.4 Relation to previous research 

Although the intrinsically infinite-dimensional character of continuous term 
structure models has been remarked by many authors, either it has been often 
dismissed as being an unrealistic working assumption in yield curve model- 
ing because of the mathematical complexity involved, leading to multifactor 
models of the term structure (Duffie & Kan [ Duffle fc Kan |, El Karoui et al. 
El Karoui et a/., 19961 |E1 Karoui fc Lacoste,1992| , |E1 Karoui et a/.,1992|l ). 



Bjork et a/. ||Bjork et a/.|| proposed a mathematical framework for arbi- 
trage pricing of interest-rate derivatives taking into account the infinite- 
dimensional nature of the yield curve process. However in their represen- 
tation, following the HJM [[HJM|| approach, this infinite-dimensional process 
is driven by a /inite-dimensional source of randomness. 



An important step was taken by Musiela's ||Musiela|| representation of the 
forward rate as a stochastic process taking its values in a space of functions. 
Musiela reformulated the HJM equations in terms of a first-order stochas- 
tic partial differential equation in the time and maturity variables. Using 
this same representation, Goldys & Musiela subsequently derived an infinite- 
dimensional version of the Black-Scholes PDE for a swaption [?]. 

Another direction which has been developed is the representation of the 
term structure as a random field with two parameters, time and maturity. 
First introduced by Kennedy [[Kennedy 1994[ , this approach allows to in- 
corporate an infinite number (or a continuum) of sources of randomness in 
the evolution of the term structure. Kennedy [[Kennedy 1994 [ considers ran- 
dom fields derived from the Brownian sheet, more general random fields are 



considered by Goldstein [[Goldstein, 1997[[ , who considers as an example ran- 
dom fields which solve a second-order stochastic PDE. Building on these 



examples, Santa Clara & Sornette [ ^anta Clara fc Sornette, 1998|| consider 



forward rate models in which the forward rate process is driven by a two- 
parameter noise process ("stochastic string shock"), again pointing out the 
relation with stochastic PDEs. 

Random field models can also be represented as multifactor models with 
an infinite number of factors. A generalization of the HJM equations to 
Gaussian infinite factor models was studied by Douady |Pouady|| . Instead of 



a random field representation, Douady uses an infinite-dimensional represen- 
tation of the yield curve and introduces a cylindrical Brownian motion as a 
source of randomness. 

An approach which unifies the infinite-dimensional character of the yield 
curve process and the possibly large number of possible sources of randomess 
was developed by Bouchaud et al. Pouchaud et al, 1997] , [Bouchaud et ai, 199^ 



based on an empirical study of deformations of the Eurodollar term structure 
between 1992 ad 1996. Bouchaud et al proposed the idea of representing 
the term structure as a randomly vibrating curve governed by a stochastic 
partial differential equation containing a second order derivative with respect 
to maturity and presented some empirical evidence in favor of such models. 

Our objective in this paper is to present a model of term structure move- 
ments which is both analytically tractable, preserves the infinite-dimensional 
character of the forward rate curve and reproduces some stylized empirical 
observations with a small number of parameters. Our approach accounts in a 
natural manner for the fact that only a small number of factors seem to gov- 
ern the covariance structure of term structure movements, without imposing 
any ad hoc unobservable state variable in the model. 

More precisely, we will try to demonstrate that specifying the term struc- 
ture process as a dynamic process with an infinite number of degrees of 
freedom and eventually proceeding to a finite dimensional approximation af- 
terwards is a more robust modeling procedure than restricting the number of 
degrees of freedom in the definition of the model, at the risk of misspecifying 
the factors. 

1.5 Outline 

The paper is structured as follows. In section ^ we recall some important 
empirical observations about term structure deformations. Based on these 
observations, we discuss in section ^ what ingredients one should incorpo- 
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rate into an interest rate model in order to reproduce the observed statistical 
properties. We then proceed to give our framework a mathematical for- 
mulation in terms of a stochastic evolution equation in a space of smooth 
functions in section ^. In the case where only local deformations are allowed, 
this equation reduces to a stochastic partial differential equation: we study a 
simple example of such an evolution equation in section and show that, albeit 
its rudimentary structure, it reproduces many properties of term structure 
deformations in a simple manner. 

2 Statistical properties of term structure de- 
formations 

As in any applied discipline, empirical observations should be the starting 
point in the construction of stochastic models in finance. Since our aim here 
is to model the dynamical behavior of the yield curve, we shall begin by de- 
scribing some important empirical facts about term structure deformations. 
The results outlined in this section we mainly refer to ||Bouchaud et ai, 1997||, 



Bouchaud et ai, 1998|| and [[Litterman fc Scheinkman]] . 



1. Smoothness in maturity: yield curves do not present highly irregular 
profiles with respect to maturity. Of course one could argue that with 
50 or 60 data points it is difficult to assess the smoothness of a curve; 
this property should be viewed more as a requirement of market op- 
erators. A "jagged" yield curve would be considered as a peculiarity 
by any market operator. This is reflected in the practice of obtaining 
implied yield curves by smoothing data points using splines. 

2. Irregularity in time: The time evolution of individual forward rates 
(with a fixed time to maturity) are very irregular. This should be 
contrasted with the regularity of forward rates with respect to time-to- 
maturity and reveals an asymmetry between the respective roles of the 
variables t and 6. 

3. Principal components: Principal component analysis of term struc- 
ture deformations indicates that at least two factors of uncertainty are 
needed to model term structure deformations. In particular, forward 
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rates of different maturities are imperfectly correlated. Empirical stud- 
Bouchaud et ai, 1997 , P^itterman fc Scheinkmaiil] uncover the in- 



les 



fluence of a level factor which corresponds to parallel shifts of the yield 
curve, a steepness factor which corresponds to opposite changes in 
short and long term rates and a curvature factor which influences the 
curvature of the yield curve. More precisely, the third principal com- 
ponent, when projected on forward rates of different maturities, shows 
a large component at maturities around one year and small coefficients 
on the two extremities of the yield curve [[Bouchaud et ai, 1997 |. 



4. Humped term structure of volatility: Forward rates of different matu- 
rities are not equally variable. Their variability, as measured for exam- 
ple by the standard deviation of their daily variations, has a humped 
shape as a function of the maturity, with a maximum at 6' ~ 1 year 



and decreases with maturity beyond one year [p3ouchaud et a/., 1997 . 
This hump is always observed to be skewed towards smaller maturities. 
Moreover, although the observation of a single hump is quite common 



[[Moraleda|| , multiple humps are never observed in the volatility term 



structure. 

5. A multivariate process: while many previous models focused exclusively 
on the short rate process, trying to represent it as a Markov process, 
recent econometric studies seem to reject the Markov hypothesis for an 
interest rate of a given maturity, the short rate in particular, pointing 
out to the interdependence between interest rates which calls for a 
multivariate approach. 

3 Modeling strategy 

What are the lessons to be drawn from these empirical observations? We will 
now try to define some criteria which a model should try to respect in order 
to give a "faithful" statistical representation of interest rate fluctuations. 

3.1 Role of the short rate 

First, the actual dynamics (as opposed to the risk- neutral dynamics) of the 
forward rate curve cannot be reduced to that of the short rate: the statistical 
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evidence points out to the necessity of taking into account more degrees 
of freedom in order to represent in an adequate fashion the comphcated 
deformations of the term structure. In particular, the imperfect correlation 
between maturities and the rich variety of term structure deformations shows 
that a one factor model is too rigid to describe yield curve dynamics. 

Furthermore, in practice the value of the short rate is either fixed or 
at least strongly influenced by an authority exterior to the market (Federal 
Reserve, central banks), through a mechanism different in nature from that 
which determines rates of higher maturities which are negotiated on the 
market. The short rate can therefore be viewed as an exogenous stochastic 
input which then gives rise to a deformation of the term structure as the 
market adjusts to its variations. It is therefore plausible from an economic 
point of view to model separately the dynamics of the short rate. 

Second, as shown by Ait Sahalia in a recent study ||Ait-!Sahaha 1997 



the short rate and the long rate (or equivalently, the short rate and the 
spread) can be reasonably described by a bivariate diffusion such as the one 
considered by ||Brennan fc Schwarz 1979|] or [ [Schaefer fc Schwartz . 



3.2 Sources of randomness 

Traditional term structure models such as [|Vasicek|, |CIK|, |HJM|| define 



implicitly or explicitly- the random motion of an infinite number of for- 
ward rates as diffusions driven by a finite number of independent Brownian 
motions. This choice may appear surprising since it introduces a lot of con- 
straints on the type of evolution one can ascribe to each point of the forward 
rate curve and greatly reduces the dimensionality i.e. the number of degrees 
of freedom of the model, such that the resulting model is not able to re- 
produce any more the complex dynamics of the term structure. Multifactor 
models are usually justified by refering to the results of principal compo- 
nent analysis of term structure fluctuations. As remarked above (Sec. |^), it 
is often observed that the first three principal components explain more than 
95% of the observed variance of forward rates, suggesting that a three factor 
model would be sufficient. However, one should note that the quantities of 
interest when dealing with the term structure of interest rates are not the 
first two moments of the forward rates but typically involve expectations of 



^For a concise review of multifactor term structure models see [Duflie & Kan| 
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non-linear functions of the forward rate curve: caps and floors are typical ex- 
amples from this point of view. Hence, although a multifactor model might 
explain the variance of the forward rate itself, the same model may not be 
able to explain correctly the variability of portfolio positions involving non- 
linear combinations of the same forward rates. In other words, a principal 
component whose associated eigenvalue is small may have a non-negligible 
effect on the fluctuations of a non-linear function of forward rates. This 
question is especially relevant when calculating quantiles and Value-at-Risk 
measures. 

In a multifactor model with k sources of randomness, one can use any 
k + 1 instruments to hedge a given risky payoff. However, this is not what 
traders do in real markets: a given interest-rate contingent payoff is hedged 
with bonds of the same maturity.^ These practices reflect the existence of 
a risk specific to instruments of a given maturity. The representation of a 
maturity-specific risk means that, in a continuous-maturity limit one must 
also allow the number of sources of randomness to grow with the number of 
maturities; otherwise one loses the localization in maturity of the source of 



randomness in the model. This point is discussed in more detail in Sec. 5.1 



3.3 The Markov property 

An important ingredient for the tractability of a model is its Markovian 
character. Non-Markov processes are difficult to simulate and even harder 
to manipulate analytically. Of course, any process can be transformed into 
a Markov process if it is imbedded into a space of sufficiently high dimen- 
sion; this amounts to injecting a sufficient number of "state variables" into 
the model. These state variables may or may not be observable quantities; 
for example one such state variable may be the short rate itself but an- 
other one could be an economic variable whose value is not deducible from 
knowledge of the forward rate curve. If the state variables are not directly 
observed, they are obtainable in principle from the observed interest rates 
by a filtering process|bl Karoui fc Lacoste,199!^. Nevertheless the presence 



of unobserved state variables makes the model more difficult to handle both 
in terms of interpretation and statistical estimation. This drawback has mo- 
tivated the development of so-called affine curve models models where one 



^Unless, of course, liquidity considerations impose the trader to do otherwise. 
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imposes that the state variables be affine functions of the observed yield curve 
Puttie fc Kan|| . While the affine hypothesis is not necessarily realistic from 
an empirical point of view, it has the property of directly relating state vari- 
ables to the observed term structure. We will try to conserve this desirable 
feature in our model. 

3.4 Continuity of term structure deformations 

Another feature of term structure movements is that, as a curve, the forward 
rate curve displays a continuous deformation: configurations of the forward 
rate curve at dates not too far from each other tend to be similar. An an- 
imation movie of the successive positions of the forward rate curve displays 
a continuous movement where the observer can follow a given point on the 
curveQ. This continuity of deformations should be properly defined in math- 
ematical terms and accounted for in a term structure model. 

3.5 Smoothness in maturity 

As already noted above, most applications require the yield curve to have 
some degree of smoothness e.g. differentiability with respect to the matu- 
rity 9. This is not only a purely mathematical requirement but is reflected 
in market practices of hedging and arbitrage on flxed income instruments. 
Market practitioners tend to hedge an interest rate risk of a given maturity 
with instruments of the same maturity or close to it. This important obser- 
vation means that the maturity is not simply a way of indexing the family 
of forward rates: market operators expect forward rates whose maturities 
are close to behave similarly. Moreover, the model should account for the 
observation that the volatility term structure displays a hump but that mul- 
tiple humps are never observed. Such an effect could be taken into account 
by the presence of a term which couples together each forward rate with its 
neighboring points on the yield curve. We will elaborate more on this aspect 
in Sec. |^. 



Redisplay Fig. 4.4 in a postscript version of this article and watch! 
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4 A stochastic evolution equation for term 
structure deformations 

Based on the above considerations, we will now proceed to describe the defor- 
mations of the term structure in mathematical form by means of a stochastic 
evolution equation, translating each of the criteria outlined above into their 
mathematical equivalents. 

4.1 Definitions and notations 

We will parametrize the evolution of the term structure of interest rates by 
the instantaneous forward rate curve (FRC), denoted by ft{0) where the 
subscript t denotes time and 9 G [6'min, 6'max] the time to maturity. Note that 



some authors (e.g. [|H JM|| ) specify the forward rate curve as a function of the 
maturity date T; our parametrisation is related to the HJM parametrization 
inEq.(|lDby 

/HjM(t,T) = f,{9 = T-t) (4) 



As remarked in ||Musiela|| , this parametrization has the advantage that the 



forward rate curve process ft will belong to the same function space (a space 
of continuous curves defined on [6'min, 6'max]) when t varies, which is not the 
case of the process /hjm(^, •) whose domain of definition [t,T] shrinks with 
time. Here 6'min is the shortest maturity available on the market and 6'max 
the longest. r(t) = ftiOmin) will be called the short rate, l{t) = ftiOram^) the 
long rate. The quantity s{t) = l{t) — r(t) is the spread. 

In most interest rate models 6'min is taken to be and 6'max = +oo but 
this is not necessarily the best choice nor even realistic. First, it is obvious 
that in empirical applications maturities have a finite span and 6'max will be 
typically 30 years or less depending on the applications considered. Second, 
the finiteness of ^max avoids some embarassing mathematical problems related 



to the ^ — > oo limit [pybvig et al., 1996|, [El Karoui et a/.,1996|| which are not 



necessarily meaningful from an economic point of view. More importantly, 
we shall see that the 6'max = +oo limit is not "innocent" : setting 6'max to 
a large but finite value can be qualitatively different from taking it to be 
infinite. 
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4.2 Decomposition of forward rate movements 

As mentioned above, given the particular nature of the short rate and the 
well known role of the short rate and the spread as two principal factors, we 
first proceed to "factor" them out of the model and parametrize the term 
structure as follows: 

Me) = rit) + sit)[Y{9)+X,ie)] (5) 

where F is a deterministic shape function defining the average profile of the 
term structure and Xt{9) an adapted process decribing the random deviations 
of the term structure from its long term average shape. With no loss of 
generality we require: 



which results in 



l^(^min) = r(^,„ax) = 1 (6) 



Xt{9^,^) = Xti9^,,) = (7) 



The process Xt{9) then describes the fluctuations of a random curve with 
fixed endpoints. We will thus call Xt the deformation of the term structure 
at time t. Factoring out the fluctuations of the first two principal components 
then means modeling separately the process {r{t), s{t)) and the deformation 
process {Xt)t>o- 

In a Gaussian framework, the uncorrelatedness of the principal compo- 
nents would entail their independence. In particular the first two principal 
components (which are roughly the spread and the short rate) would be in- 
dependent from the deformation process Xt. We will use this assumption as 
a working hypothesis: 

Assumption. The deformation process Xt is independent from the short 
rate r(t) and the spread s(t). 

4.3 The short rate and the spread: a bivariate Markov 
process 

As in IBrennan & Schwarz 197911 , one can consider that the short rate and the 



long rate (or equivalently, the short rate and the spread, see [^chaefer fc Schwartz]] ) 
are well described by a bivariate diffusion process: 
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drt = fXi{rt,St)dt + ai^i{rt,St)dWl + ai^2{rt,St)dW^ (8) 

dst = H2irt,St)dt + a2Airt,St)dWl + a2,2{rt,St)dWt (9) 

where W^,W'^ are two independent Wiener processes. This formulation is 
empirically motivated by the econometric studies refered to above [Ait-Sahalia 1997| 



in which it was shown that the hypothesis of a bivariate diffusion for (r^, st) 
is not rejected by non-parametric tests while it is rejected for rt taken indi- 
vidually. 

However, the only hypothesis we need here is the jointly Markovian char- 
acter of {rt,St); ,for example, the noise source in Eq.(|^) could be replaced 
with a non-Gaussian Levy process without modifying what follows. 

4.4 Term structure deformations as Markovian curves 

We are now left with the deformation process {Xt)t>o to model. The first 
requirement we impose on Xt is its smoothness in maturity: at a given 
time t, Xt is a function defined on [6'min, 6'max] determined by the forward 
term structure which, as remarked above, is a "smooth" function of the 
time to maturity 6. Xt should therefore belong to a suitable space H of 
smooth functions which will then be the state space of our model. In view of 
interpreting our results in terms of principal component analysis, we would 
like the state space H to have some Hilbert-like structure in order to define 
orthogonal projections of Xt onto a suitable basis of H. 

The second requirement we impose is that Xt be a Markov process in 
H. This property, as remarked in |[Musiela|| , is already verified in the Heath- 



Jarrow- Morton |[HJM|| framework for the forward curve process ft- Here 
we require slightly more, namely that the Markovian character respect the 
factorial decomposition (|^). That is, we require the endpoints (r(t), s(t)) and 
the deformation Xt to be separately Markovian. 

In a functional space, there are of course a wide variety of Markov pro- 
cesses. It is the hypothesis of continuity in time of the deformation process 
(see ( |3.4| ) which enables to single out, the only class among all Markov pro- 
cesses having this property, namely diffusions. More precisely, stating that 
Xt is a if -valued diffusion process means that there exist a drift functional b 
and a volatility functional cr, defined on H, such that the evolution of Xt is 
given by a stochastic differential equation in H (written here in Ito notation): 
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dXt = fx{Xt)dt + aiXt)dBt (10) 

where Bt is an appropriate generalization of Brownian motion taking values 
in H. Here /i and a are allowed to depend on the contemporaneous term 
structure i.e it can be a function of the whole curve Xt{6), 9 G [0, 9*]. 

Formally, Eq.(|IO|) is a stochastic differential equation in an (infinite- 
dimensional) functional space H. In order to give a proper meaning to 
Eq. (p!0|) , one should start by specifying the nature of the random noise source 
Bt such that the stochastic integral implicit in Eq.(p!oD can be properly de- 
fined. There are several ways to define a generalization of the Wiener pro- 
cess and a stochastic integral in an infinite dimensional space. The relation 
between these different constructions was clarified by Yor [ Yor 1974 who 



showed that the natural setting for constructing infinite dimensional diffu- 
sions is a Hilbert space. Given a (separable) Hilbert space H (for exam- 
ple H = L'^{[9rain,0max\,^) foi somc mcasurc z/), one can define a cylindri- 
cal Brownian motion on i^ as a family {Bt)t>o of random linear functionals 
Bt : H ^ R satisfying: 

1. y(j)eH,Bo{(j)) = 

2. \/(f> E H, Bt{4>) is an J^t - adapted scalar stochastic process. 

3. Wcj) E H — {0}, -4^ is a one-dimensional Brownian motion. 

In particular, if one takes any orthonormal basis (e„) in H then its image 
{Bt{en))t>Q form a sequence of independent standard Wiener processes in R. 
This property is useful for building finite-dimensional approximations. 

A suitable choice of state space verifying such requirements is a Sobolev 
space H"^, namely the space of functions g G L'^{[9mm,9ma.x],^) such that 
the s-th derivative g^^"* is also in L'^{[9mm,9ma.x\,i^), for some measure u on 
[^min? ^max]- Givcu that the derivative df/d9 is assumed to exist in many 
applications, we would like to require s > 1. See also the discussion of this 
point in [ pouadyj ] , where a similar choice is adopted in a slightly different 
framework. 
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Figure 1: Various configurations adopted by the function Xt{9) defined in 
Eq. 1^ as obtained from Eurodollar future contracts (1992-7). This figure 
may be seen as a visualization of the stationary density of the process Xt. 
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5 Local deformations, stochastic PDEs and 
string models 

At this level of generality, not much can be said of the properties of the 
solutions of Eg . ([To|) . In this section we will show how the description of term 
structure deformations through level, steepness and curvature of the yield 
curve reduces Eq. (|lOD to a stochastic partial differential equation, of which 
some simple examples are given. 

5.1 Market segmentation and local deformations 

As mentioned before, the maturity 9 is not simply a way to index different 
forward rates and instruments: the fact that fixed income instruments are 
ordered by maturity is important for market operators. For example, this is 
reflected in the hedging strategies of operators on the fixed income market: to 
hedge an interest rate risk of maturity ^ = 8 months, an operator will tend to 
use bonds (or other fixed-income instruments) of maturity close to 8 months: 
6 months, 9 months. Although this strategy seems quite sensible, it does not 
correspond to the picture given by multifactor models: in a k-factor model, 
any k + 1 instruments can be used to hedge an interest rate contingent claim. 
For example in a two factor model one could use in principle a 30 year bond, 
a 10 year bond and a 6 year bond to hedge an instrument with maturity of 
two years! Needless to say, no sensible trader would follow such a strategy, 
which shows that in practice the factors which explain 95% of the variance 
are not enough to hedge 95% of the risk of an instrument with a non-linear 
payoff: this is precisely our principal motivations for introducing maturity- 
specific sources of randomness i.e. one independent source of randomess per 
maturity. 

The existence of maturity-specific risk naturally leads to a market for 
such risk. Indeed, some macroeconomic theories of interest rates have con- 
sidered the interest rate market as being segmented : for example in a first 
approximation one can consider the market for US Treasury bills. Treasury 
notes and Treasury bonds as being 3 separate markets where prices are fixed 
independently. In a continuous-maturity model this would mean that the 
interest rate market is partitioned into independently evolving markets in- 
volving instruments with maturity between 9 and 9 + d9. 
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However this is not strictly true: as shown by principal component anal- 
ysis, long-term rates react to variations in the short rate in a way that is 
not explainable simply via parallel shifts and vertical dilations of the term 
structure. One way to conciliate the interdependence of rates of various ma- 
turities with the segmentation of markets across maturities is by considering 
deformations of the term structure that are local in maturity: a forward rate 
of maturity 6 is more sensitive to variations of rates with maturity close to 
6. We are not dealing here with a strict segmentation of the market into 
separately evolving markets but a "soft" segmentation which simply implies 
that the market for each maturity adjusts itself to the variation in rates of 
maturities immediately above and below it. This means for example that, 
among all rates of maturity > 1 year, the 1 year rate will have a higher 
sensitivity and react more quickly to a variation in the short rate since it is 
closer in maturity. 

5.2 Level, steepness and curvature 

In mathematical terms, the local deformation hypothesis means that the vari- 
ation of Xt{6) will only depend on the behavior of Xt{.) around 6. How can 
one parametrize the shape of the term structure around a given maturity 67 
Given that Xf : 6 ^ Xt{6) is assumed to be a smooth curve, its local behavior 
around 6* can be described by its first few derivatives: Xt(^), deX^, dgX^, 

As noted before, empirical studies seem to identify the level of interest 
rates, the steepness (slope) of the term structure and its curvature as three 
significant parameters in the geometry of the yield curve [[Litterman fc Scheinkmaii 



In a market involving instruments of maturity between 9 and 9 + d9, these 
three features are described by the level of rates, and the first two derivatives 
with respect to 9. Combining the local deformation hypothesis formulated 
in Sec.^TT] with a local description of the term structure by level, steepness 
and curvature one obtains that the drift and volatility of Xt{9) can only de- 
pend on Xt{9), dgXt, d^Xt. Therefore, Eq.(|ID|) becomes a second order 
stochastic partial differential equation: 

dx, ^ [^ + Kxm. '-§. ^)] dt + aixm, '-§. '^)dBm 

Vt > 0,X,(^^,J = Xt{9^^) = Xi=o(^) = ^o(^) (11) 
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This equation is the mathematical expression of the fact that deformations 
are local in maturity and that the deformation at maturity 9 depends on the 
level, steepness and curvature of the term structure around 6. 

In the general case where h and a are smooth but nonlinear functions 
of their arguments, Eq.(|ll]) is not easy to study: indeed, it is not trivial 
to define properly what is meant by a solution of Eq.([TT|) and even less to 
study their regularity. An approach to the fully non-linear case using the 
notion of viscosity solution has been recently proposed for the case of a 
noise-source depending only on t [[Lions fc Souganidijl . In order to point out 



the differences with HJM-type models resulting from the local deformation 
hypothesus, we shall consider the case of a forward rate dynamics such as 
(|T]) which is perturbated by a term depending on the curvature: 

dX, = [-g^ + hit, e, Xie)) + 2 ^^] ^^ + ^(^' ^' X^{Q))dB,{Q) 

Vt > 0,Xi(e^i„) = Xi(e^ax) = X,=o(^) = ^o(^) (12) 

Properties of stochastic PDEs such as Eq. (|1^) have been studied by Pardoux 
[Fardoux| and collaborators. Properties of term structure deformations de- 
scribed by Eq.([T^) are studied in [|Cont 1998|| . In the following section we 
study the simplest case where volatility is constant; surprizingly, we will 
show that this simple case already presents many of the desirable features 
enumerated in Sec.|3|. 



5.3 The linear parabolic case 

In order to illustrate what are the type of dynamics implied by Eq . (|T^) 
for term structure deformations, we will now study the simplest example of 
the above equations which incorporates the influence of local steepness and 
curvature, namely the case where a is independent of X^. For the sake of 
simplicity we will deal here with the constant volatility case but all the results 
below remain valid in the case of an arbitrary deterministic function of time 
t (for details see ||Cont 199^ ). The case of constant volatihty leads us to the 



following stochastic partial differential equation 



5 



dX ,dX Kd^X^, ,„,^, 



^Up to the first derivative in 6, this closely resembles what is known as the stochastic 
heat equation in the PDE literature. 
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Vt > 0, XtiO^i,,) = Xi(^„,ax) = (14) 

V^ e [^^in, ^„.ax], Xt=oiO) = Xoi9) (15) 

5.4 Eigenmodes and principal components 

Let 9* = ^max — 6'min be the maturity span of the observed forward rate curve. 
By translating the maturity variable one can assume ^mm = without loss 
of generality in what follows. We consider as state space for our solutions 
the Hilbert space H of real-valued functions defined on [0, 6*] with the scalar 
product: 



/■"* 2/9 

<f,g> = d9 exp(-)/(%(^) (16) 

Jo K, 

The subscript H in < ., . >h will be omitted in most of this section. Let A 
be the operator in H defined by: 

du Kd'^u ,,„, 

It is not difficult to show that A has a discrete spectrum, with eigenvalues 
and eigenfunctions given by: 

A.Cn = -A„e„ (18) 

/ 2 ,n9n^ , 9 ^ , , 

^n{0) = ^j0;M^)eM--) (20) 

where n takes all integer values > 1. Here the eigenfunctions e„ have been 
normalized such that (e„)„>i is an orthonormal basis of H. 

The functions (eigenmodes) e„(^) play the role of the principal components 
for the deformation process Xt. That is, if we perform a principal compo- 
nent analysis on a realization of the process Xt, for a large enough sample the 
empirical principal components would reproduce the eigenmodes en{9). The 
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first two of these eigenmodes are sliown in Fig.( |5.4| ). Tlie role of tlie exponen- 
tial term in Eq.(pOD is clearly visible: the eigenfunctions become "skewed" 
towards shorter maturities and only a single hump, whose position is deter- 
mined by the value of k, is visible. Recall that this exponential term stems 
simply from the fact that we are parametrizing the forward rate process by 
time to maturity 9 instead of maturity date T [[Musiela]| . In particular, in 



contrast with multifactor models ||Moraleda |, there is no need to use a com- 



plicated volatility structure o"(t, 9) to obtain a volatility hump. The position 
of the hump gives a (first) simple method for calibrating the value of k to 
empirical observations. 

The Green function (propagator) associated to the operator A may then 
be expressed in terms of an eigenmode expansion: 

G(t,x,y) = ^ exp(-A„t)en(a;)e„(y) (22) 

n>l 

and Eq.(|T3p can be properly defined in the following integral form: 

Xt{9) = f G{t,9,y)Xo{y)dy+ f ds f G{t - s,9,y)a,dB,{y) (23) 

JO JQ Jo 



Let Xt{9) be the solution of Eg. (|23D . Define the coordinates of the solution 
in the eigenvector basis as: 

r^* 29 

Xnit)=<Xt,en> = / d9 Xt{9)e^{9) exp(— ) (24) 

Jo K 

The coefficients x„(t) therefore represents the projection of the deformation 
process Xj on the n-th principal component (eigenmode). For each n, x„(t) 
is then a solution of a linear stochastic differential equation: 

dxnit) = -XnXnit)dt + aodWl' (25) 

where IV" = Bt{en) are independent standard Wiener processes. The x„(t) 
therefore consitute a sequence of independent Ornstein-Uhlenbeck processes: 

x^{t) = e-^"' <en,Xo>+ f ds e'^-'^'-'^aodW^ (26) 

Jo 

This last equation has an interesting interpretation. Remember that a;„(t) the 
projection of the deformation process Xf on the n-th principal component. 
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Eq.(|26|) expresses x„(t) as the sum of two components, the first one being 
the contribution of the initial term structure to Xn and the second one its 
stationary value. Eq.(^B|) may then be interpreted by stating that the forward 
rate curve "forgets" the contribution of the n-th principal component to the 
initial term structure at an exponential rate with characteristic time 



A„ 1 + 



KZT- 



Therefore, a perturbation of the initial term structure due to the n-th prin- 
cipal component will disappear or be smoothed out after a typical time r„ 
which decreases with n: "singular" perturbations die out more quickly than 
smoother ones. This shows the important relation between the smoothing 
property of the deformation operator A and the decay of its eigenvalues: 
an operator with a quickly decaying spectrum will guarantee a fast decay 
(in time) of the singularities appearing in the term structure and restore 
smoothness in maturity. This gives a second interpretation of the parameter 
k: in addition to determining the position of the volatility hump, it also 
determines the decay rate of perturbations of the term structure. This in- 
terpretation gives a second, independent method for calibrating the model 
parameters to empirical data. One can easily imagine more general models 
where these two roles of n can be attributed to two separately calibrated 
parameters [|Cont 1998 |. 



The Ornstein-Uhlenbeck process is the process used to represent the short 



rate in |[Vasicek|| : it possesses the fundamental property of mean reversion 
which has made it a popular model in interest-rate modeling. In our case, 
this mean reversion is observed in the principal components of the yield 
curve: while individual forward rates may have a non-stationary and irregular 
behavior the yield curve as a whole will converge to a stationary state with 
a mean-reverting behavior. 

Given the explicit form of A„, it is easy to show that the series 



E[Xt{e)] = Y.E[xn{t)]en{e) (28) 

n>l 

Var[Xt{e)] = ^rar[x„(t)]e„(^)2 (29) 



n>l 
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First eigenmode 
Second eigenmode 
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Figure 2: First two eigenmodes of the operator A, with K=4 and 9* = 30 
years. Note the maxima situated at short maturities. 



are absolutely convergent for all (t, 6) and the sum 



(30) 



n>l 



defines a unique Gaussian random field with mean and variance given by: 



E[Xt 
Var[Xt 



n>l 



T^Ant 



< Xo,e„ > e„(6') 



CTr 



;E 



1 — exp(— 2A„t) 

2 An 



en{ey 



(31) 
(32) 



5.5 Average term structure and mean reversion 

Under the above assumptions, one can calculate the average shape of the 
term structure of forward rates from (Eq. ^). 



E[f{tM= EW)]+ E[s{t)]Y{d) 



(33) 
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The shape function Y{6) can therefore be chosen in order to reproduce 
the average term structure. In jBouchaud et al, 1997|| it was found that the 



function 



yi(^) = \hr (34) 



gives a good fit of the average shape of the Eurodollar term structure for 
maturities ranging from 3 months to 10 years. However, the precise analytic 
form of the shape function Y does not affect the results above. How does 
the yield curve fluctuate around its average shape? It is easily seen from 
Eq. (|3l|) that the process Xt converges to a Gaussian random field X^ with 
mean zero and covariance: 

Covix,ie),xAe')) = ^lY^^JB^ji^Kf^^l (35) 

n>l ^^" 

In terms of term structure movements, stationarity of term structure de- 
formations implies a mean-reverting behavior of forward rates. But a model 
such as the above asserts more: it enables to calculate (therefore calibrate, 
if one is interested in such quantities) the probability of a given yield curve 
deformation. 

5.6 Random strings: parabolic vs. hyperbolic formu- 
lation 



In a recent work [panta Clara fc Sornette, 1998|| it has been proposed to con 



sider stochastic partial differential equations of hyperbolic type to describe the 
evolution of the forward rate curve. Such an equation differs from the above 
one through the presence of a second-order time derivative which dominates 
the dynamics. The model in ||Santa Clara fc Sornette, 1998|| is formulated 



in terms of the forward rate process itself (in the spirit of | [HJM| |) and not 
in terms of the deformation process Xj, by examining the effect of inserting 
a second-order time derivative in the equations of Sec. |[ Let us there- 
fore consider the general case where the evolution equation contains both a 
propagation term and a diffusion term: 

B^x dx dx Kd^x ,„,^, 
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Vt > 0, Xi(0) = Xt{e*) = (37) 

r)X 

x,=oW = Xo{e) — ^^ = yo (38) 

The above equation is analogous to that of a vibrating elastic string, hence 
the name of "string models" given to such descriptions of term structure 
movements. The case / = is the one studied in Sec.|^; the case f ^ oo 
(the other parameters being appropriately rescaled) is the stochastic wave 
equation ||Da Prato fc Zabczyk|| , the "space" variable being 6 + t. Note that 



the stochastic PDE used in [^anta Clara fc Sornette, 1998|| is formulated as a 



PDE perturbated by a two-parameter (" space-time" ) noise (called a "stochas- 
tic string shock") while our Eq.(|12|) or the general case Eq.([Tl|) was presented 
above as an evolution equation for a curve in some function space. In the 
case of a parabolic SPDE, where only the first derivative with respect to time 
is involved, the two approaches are equivalent for a two-parameter process. 
Choosing one approach or the other then amounts to viewing the solution 
of a stochastic PDE either as a random field or as a stochastic process in 
a function space. Implicit in this choice is whether the object of interest 
for modeling purposes is an individual interest rate or the deformation of a 
multivariate object, namely the term structure. The parabolic equation in 



Sec. |5.3| has the merit of emphasizing the asymmetric roles of the variables 9 
and t, an empirically desirable feature which is not present, as we shall see 
below, in [ [Santa Clara fc Sornette, 1998 |. 



Manifestly, the operator on the right hand side of Eq.(|36D is the same as 
in Eq.(|T^). This means that the deformation eigenmodes (the eigenf unctions 
of A) will remain the same as in Eq.([T3|) studied above but the projection of 
the process Xt on each of them will be different. For example, a stationary 
solution of Eq. ( |36[ ) will not give the same weight to the eigenmodes as and 
therefore the results of a Principal Component Analysis of Eq . (|36D will differ 
from that of Eq.(|r^) in terms of the eigenvalues. 

The representation of the equation in terms of its projections on the 
eigenmode basis e^ gives as above a stochastic equation for the scalar process 
Xn(t): 

f^ + ^ = -Kyn{t) + aoWr (39) 

This formal second-order stochastic differential equation is interpreted in the 
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usual way, as follows. Consider the Green function Unit) of the operator 
9*2 "•" at 



f^ + lW + ^nu, given by: 



U4t) - ^-,=_[e---e-1l.>o (40) 

(41) 

where ri „ and r2^n being the roots of the associated characteristic equation: 

fr^ + r + Xn = 0. (42) 

The process x„(t) is then given by the stochastic convolution integral: 

yn{t) = ane-'-^-' + bne-'-'-' + JUn{t-s)dW: (43) 

which is well-defined since the integrator is square- integrable in s. Here a^ 
and bn are defined by the initial term structure. Depending on the values of 
/ and K,, two scenarios are possible: 

1. Oscillating initial conditions: if 

/ > ^^^ (44) 

'2(1 + ^) ^ ^ 

then for all n > 1, Eq.^ has two complex conjugate roots given by: 

zVAXJ- -1-1 1 . 

ri,n = 2/ " "2/ "^'^ ^^ 

-tV4Kf -1-1 1 . 

r2,n = 2/ " "2/ ~ ^^" ^ ^ 

The real part gives an exponential damping of the initial conditions 
which characterizes the mean reverting behavior of Xj as in the parabolic 
case. First remark that, unlike the parabolic case where "bumpy" prin- 
cipal components which contribute the most to non-smoothness in ma- 
turity decay more quickly, here all principal components decay with 
the same speed i.e. a mean reversion time of 2/. Recall that n still 
determines the position of the volatility hump so k ~ 1 year. So (^) 
implies that f > Q months. The mean reversion time of the whole 
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curve is thus around a year. However, a new phenomenon appears: the 
principal components do not simply revert to their mean but oscillate 
around their mean with a frequency C(j„/27r which increases with n: 

Xnit) = A„e-*/2^ cos(cu„t + 0„) (47) 

The phase 0„ and amplitude factor A„ are determined by (two) initial 
conditions (see below) . The oscillation of the term structure around its 
mean is not necessarily an undesirable feature of this model and indeed 
can be justified on economic grounds fCont 1998| . But the slow mean 



reversion combined with increasingly faster oscillations of the higher 
order principal components leads to non-smoothness in maturity of the 
solutions of Eq.(^) |Pa Prato fc Zabczyk|| : in fact one should expect 



the cross-sections in time or maturity to have the same irregularity 
which, as pointed out in Sec.^ is not a desirable feature for a term 
structure model. 

2. Selective damping of principal components: if 

/ < ^^^^ (48) 

2(1 + ^) ^ ^ 



then 3A^ > 1 such that for n < N Eq . (^2]) has two real, negative roots 
whereas for n > N the roots are complex conjugates with negative 
real parts. The projections of the deformations process Xt on the first 
A^ eigenmodes will have a mean reverting behavior as in the parabolic 
caseQ, with a mean reversion time increasing with n. For n > N, 
Xn{t) will have a damped oscillatory behavior, with a damping time 
r = 2/ independent of n and an oscillation frequency increasing with 
n as above. 

Another crucial difference between Eq.(^) and Eq.([T^) is the nature of 
the initial conditions. In the case of the parabolic equation (0) the problem 
has a well defined solution once the initial term structure is specified through 
Xq. This is not sufficient in the case of Eq.(^): one must also specify the 
derivative with respect to time at t = 0. In the case of a vibrating string, this 



^In fact the decay of the initial condition is described in this case by the superposition 
of two decreasing exponentials with time constants given by rj^ and r^„ 



30 



means specifying the initial position and the initial velocity of each point of 
the string. For a model of the forward rate curve, this can be inconvenient: 
while the initial term structure is the natural input for the initial condition 
of a dynamic model, the time derivative of the forward rates is not easily 
evaluated, especially given the irregularity in time of forward rate trajectories 
which prevents such a model from being calibrated in a numerically stable 
manner 



Cont 1998 



We therefore conclude that the question of including a second-order time 
derivative in Eq. (|36D is not simply a matter of taste: the presence of a second 
derivative radically changes both the dynamic properties of the equation and 
the nature of the initial conditions needed to calibrate the model, in an 
empirically undesirable fashion. Our analysis thus pleads for a description of 
yield curve deformations through a parabolic rather than hyperbolic SPDE. 

6 Conclusion and perspectives 

We have presented a simple stochastic model for describing the fluctuations 
of the term structure of forward rates: the forward rate curve is described as 
a random curve oscillating around its long term average. The model studied 
in Sec. ^ should be viewed as the simplest example of the type of model pre- 
sented in Sec. |. However this simple example has the benefit of emphasizing 
the role of the second derivative with respect to maturity in the evolution of 
the term structure: indeed, as we have seen above, it is this second deriva- 
tive which tames the potentially infinite number of sources of randomness 
and maintains a regularity in 9 while allowing for independent shocks along 
maturities. It also gives the correct form for the principal components as well 
as a qualitatively correct estimate for their associated eigenvalues. These re- 
sults show the importance of the concept of local deformation explained in 
Sec.|0|, of which our equation is the simplest example. The model in Sec.p75| 
can be easily generalized to the case where the volatility surface a{t, 6) is 
an arbitrary deterministic function [JUont 19l?E|. The introduction of a non- 



linear drift h depending on the level of interest rates - as in Eq. ([T2|) - is also 
possible. More general cases remain to be studied. 

As mentioned in the introduction, our objective has been to obtain a faith- 
ful continuous-time representation of the statistical properties of the forward 
rate curve. What remains is to establish the link with the arbitrage pricing 
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approach and examine the constraints imposed by absence of arbitrage on 
models of the type exhibited above. Previous work in this direction [ Douady 



indicates that such an analysis requires a careful reconsideration of the class 
of arbitrage strategies one is willing to consider in the context of a partially 
segmented market such as the fixed- income market. 
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